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ABSTRACT 

We investigate the eorrelation of star formation quenching with internal galaxy proper- 
ties and large scale environment (halo mass) in empirical data and theoretical models. 
We make use of the halo-based Group Catalog of Yang and collaborators, which is 
based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer 
(GALEX) are also used to extract the recent star formation rate. In order to investi- 
gate the environmental effects, we examine the properties of "central" and "satellite" 
galaxies separately. For central galaxies, we are unable to conclude whether star forma- 
tion quenching is primarily connected with halo mass or stellar mass, because these two 
quantities are themselves strongly correlated. For satellite galaxies, a nearly equally 
strong dependence on halo mass and stellar mass is seen. We make the same compar- 
ison for five different semi-analytic models based on three independently developed 
codes. We find that the models with AGN feedback reproduce reasonably well the 
dependence of the fraction of central red and passive galaxies on halo mass and stellar 
mass. However, for satellite galaxies, the same models badly overproduce the fraction 
of red/passive galaxies and do not reproduce the empirical trends with stellar mass or 
halo mass. This satellite overquenching problem is caused by the too-rapid stripping 
of the satellites' hot gas halos, which leads to rapid strangulation of star formation. 



1 INTRODUCTION 



Galaxies may be broadly divided into two categories: those 
that are forming stars fairly rapidly relative to their past av- 
eraged star formation rate (SFR), and those that show little 
recent star formation relative to their past average. It is 
well known that, at least at low redshift {z ^ 1), the former 
type tend to have blue colours and to be morphologically 
disk dominated, while the latter tend to have red colours 
and to be morphologically early type or spheroid domi- 
nated. Galaxies exhibit col our bimodaliti cs thro ughout a 
wide range of cosmic history ([Strateva et aI.ii2001.:,Bell et al 



2004 iBaldrv et al.ll2004l : iBalogh et al.ll2004l : iBlanton et al 



more closely correlated with internal galaxy properties such 
as mass or luminosity, or large scale environment (sometimes 
referred to as "Nature or Nurture")? 

There have been many studies on the impact of envi- 
ronme nt on galaxy properties. For example, iDavis fc Gelled 
and many others found that early-type galax- 
ies are more s trongly clustered than the late- types, and 
iDressleJ (|198G| ) systematically demonstrated that the frac- 
tion of elliptical and SO galaxies is higher in denser en- 
vironments. Similarly, it is well known that galaxies in 
dense environments te nd to be red and to have depressed 
star formation rates (iHashimoto et aL Lewis et al 



2005b|). One of the fundamental questions in galaxy forma- 
tion is: what are the main physical forces that regulate and 
in some cases quench star formation? Are these processes 



200: 



Gomez et al]|2003l: iTanaka et al.! 20041: iBalogh et al 



2004: iKauffmann et all |2004| : iGhristlein fc Zabludofj |2005| : 
Poggianti et al.l 20061 ). There are several physical pro- 
cesses associated with environment that may play a 
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role in galaxy transformation. Major (near equal-mass) 
mergers can transform spi ral galaxies into ellipticals 
ijToomre fc Toomrd [l973 : cf. iBarnesI [200I 1 . and may also 
quench future star formation by ejecting the ISM from 
the galaxy via s tarburst, AGN, or shock-driven winds 
dCox et all I2OO4I: ISprineel. Pi Matteo. fc Hernauisd I2OO5I : 
[Murray. Quataert. fc ThompsonI l2005h . In rich clusters , 
where the probability of merging is suppressed because of 
the large relative velocities of galaxies, galaxy "harrassment" 
(rapid encounters or fly-bys) may cause a less dramatic 
form of transformation by heating disks, perhaps causing 
the for mation of a bar and growt h of a spheroidal com- 
ponent iMoore. Lake, fc Katj 0*9981 ). Also, cold gas can be 
stripped out of the galaxy both by tidal forces due to the 
background dark matter dominated potential of the clus- 
ter, and due to ram pressure stripping b y the intracluster 
medium (iGunn fc Gott 1972: Abadi. Moore, fc BoweJI 19991 : 
iQuilis. Moore, fc Boweil[2000l : IChung et aLlbOOTl ). Sirnilarlv. 
the hot halo that provides future fuel for cooling and star 
formation may be efficiently stripped in dense environ- 
ments, thus quenching f urther star formation by "starva- 
tion" or "strangulation" jLarson. Tinslev. fc Caldwelll Il980l : 
iBalogh fc MorrislbOOdrBekki. Couch, fc Shioval l2002l) 

However J studi es by Ivan den Bosch et all (1200881 ) and 
iTanaka et all (|2004 ) suggest that processes specific to clus- 
ters (e.g. ram-pressure stripping) are not the main mecha- 
nisms for quenching star-formation act ivity. Similar result s 
were also found at higher redshift (e.g.. lCooper et al]|2006h . 
Moreover, both the morphological and spectrophotometric 
characteristics of galaxies are also known to be strongly 
correlated with their internal pro perties, such as luminos- 
ity, mass, and interna l velocity (jRoberts fc HavnesI 1 19941 : 
iKauffmann et al.ll2003l ). More massive or luminous galax- 
ies are more likely to be spheroid dominated, to be red, and 
to have old s tellar populations and li ttle recent star forma- 
tion. Indeed, IKauffmann et al] (|2003l ) showed that galaxies 
appear to make a transition in all of these properties above 
a critical stellar mass of ^ 3 X 10^° A/0. 

Again, there are various physical processes that one 
might expect to imprint this kind of dependence on inter- 
nal properties. Supernova feedback has long been invoked 
as a mechanism t hat could heat and drive gas out of galax- 
ies (|Larsonlll973 : lDekel fc Silklll986l ). and is expected to be 
more effective in low-mass galaxies. There is also mounting 
observational evidence that AGN are associated with the 
quenching of star formation (|Schawinski et al]l2006l . l2007bl : 
ISalim et al.ll2007l ). AGN feedback is expected to have more 
impact on massive galaxies , which host larger mass black 
holes (e.g. lSilk fc Re3ll998l ). 

The emerging picture is that AGN seem to have 
two modes of fueling and also to couple to their sur- 
roundings in different ways. "Bright mode" AGN are 
associated with high (near Eddington) fueling rates, and 
observationally with classical X-ray or optically bright 
quasars. The radiation emitted by these objects can couple 
with the cold gas in the galaxy, perhaps driving power- 
ful winds that can drive the gas out of the galaxy and 
quen ch sta r formation jMurrav. Quataert. fc ThompsonI 



20051: 



Di Matteo. Springel. fc HernquisT l2005l : 



Monaco fc Fontanotll2005l ) 



ton rate, and which typically do not show classical quasar- 
like signatures such as bright X-ray radiation or broad emis- 
sion lines. However, these objects are associated with the 
efficient production of radio jets, which may be able to 
couple to and heat the ho t gas in the gala:xy's halo (e.g. 



iMcNamara fc Nulsen|[2007 ). These objects are often referred 



to as "radio mode" AGN ( Croton et al.|[200^ ). This process 



is expected to be a function both of internal properties and 
environment: the bigger the black hole, the more energy can 
be tapped (and black hole mass is of course correlated with 
galaxy mass) , and empirically it is known that the fraction of 
radio-detected ga laxies increases strongly with stellar mass 
(|Best et al.ll2005l ) and halo mass (|PasauaU et al.ll2008l ). But 
as well, the radio jets must have a "working surface" and 
therefore can only be effective in environments that can sup- 
port a quasi-hydrostatic hot gas halo, such as groups and 
clusters. Galaxies in smaller mass halos (Mhaio ;S 10^2 Mq) 
likely accrete m ost of their gas in a "co l d ffow", and never 
form a hot halo dBirnboim fc Dekelll2003l :[ Dekel fc BirnboimI 
l2006l : lKeres et al.ll2005l ). 

A great deal of recent progress has been made 
towards developing a comprehensive theory of galaxy 
formation. The semi-analytic approach, although it has 
its limitations, is a powerful and flexible tool for ex- 



White fcFrenkI 1991: 'Kauffmann. White, fc Guiderdonil 


19931: IKauffmann et all Il999: SomerviUe fc Primackl 


I999I; 


Cole et all I2OOOI: SomerviUe. Primack. fc Fabeij 


200 1|: 


Springel et al.l I2OOI: iBenson et al. 


20031: Hatton et al.1 


20031: Khochfar fc Burkerti |2003|. 12005 


: Kang et al. 20051)? 



Several groups have now implemented one or both modes 
of black hole growth an d AGN feedback in t o their semi- 
analytic models (e.g.. [ Croton et al.l l2006l: iBower et al 



am 

2009: iMonaco. Fontanot. fc Taffonil l2007l : ISomerville et al 
2008bl ). There seems to be consensus that including these 
new processes leads to greatly improved agreement with key 
observations such as galaxy luminosity or mass functions, 
and the galaxy colour-magnitude distribution or stellar 
mass vs. specific star formation distribution. However, this 
necessitates including several new recipes and parameters 
associated with the poorly understood physics of black 
hole growth and AGN feedback. Each model contains 
somewhat different parameterizations and treatments of 
these processes, yet they all produce similar results for 
global quantities such as the galaxy luminosity function and 
colour-magnitude distribution, no doubt in part because 
these observations were "targets" that the modellers were 
trying to reproduce. One goal of our work here is to 
determine whether breaking down the fraction of quenched 
galaxies in the dual space of internal galaxy properties 
and environment can discriminate between these different 
treatments of AGN feedback. 

The large and homegeneous databases provided by 
modern surveys, such as SDSS and GALEX, finally al- 
low a statistically significant investigation of different sub- 
populations of galaxies. However, an obvious question that 
arises in any study of galaxy environment is exactly how to 
measure and characterize environment in observational sam- 
ples. Clearly, it is desirable to span as broad a range of envi- 
ronments as possible, from isolated field galaxies to groups 
to rich clusters. The majority of the studies in the literature 
have parameterised environment in terms of the number of 



In contrast, many massive galaxies seem to contain 
AGN which are accreting at a small fraction of their Edding- 
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galcixies, either within a fixed metric aperture, or by the 
distance of the rith-nearest galaxy, where n is typically in 
the range 3-10. Although these indicators are straightfor- 
ward to measure, they are not straightforward to interpret 
in physical ter ms or to compare with t heor etical models (see 
discus sions in iKauffmann et al] |2004| and IWeinmann et all 
l2006al ). An alternate approach is to use a galaxy group cat- 
alogue, in which galaxies in an observational catalog are not 
only grouped together into putatively gravitationally bound 
structures, but also the total mass of their associated dark 
matter halo is estimated (Yang et al. 2007). Thus, the rela- 
tionship between galaxy properties and dark matter halo 
properties can be studied directly. Another advantage of 
this method is that galaxies can be separated into "cen- 
tral" and "satellite" populations. Most of the environment- 
related tranformation mechanisms described above (such as 
stripping) are expected to work only on satellite galaxies, 
so this offers a way to separate out the effects of different 
physical processes. 

Here we make use of a large galaxy group catalogue con- 
structed from the S PSS using th e halo- based galaxy group 
finder developed by lYang et al.l (|2005h . These catalogues 
have already been used for several studies r egarding the envi- 
ronmen t dependence of galaxy properties. IWeinmann et all 
(|2006al ). using the version based on the SDSS DR2, studied 
the correlations between colours, specific star formation rate 
and halo mass. Splitting the galaxy population into early 
and late types, based on their colours and specific star for- 
mation rates, they found that, at a fixed luminosity, the late 
(early) type fraction of galaxies increases (decreases) with 
decreasing halo mass. Using the much larger galaxy group 
catalogue of Ya ng et al. (2007; hereafter Y 07), based on 
the SDSS DR4, Ivan den Bosch et al.l (|2008ah showed that, 
on average, satellite galaxies are redder and more concen- 
trated than central galaxies of the same stellar mass. They 
also found that the colour and concentration differences of 
central-satellite pairs matched in stellar mass are completely 
independent of the mass of the host halo of the satellite 
galaxy. This indicates that satellite-specific transformation 
mechanisms are equally efficient in host haloes of all masses 
and rules against satellite transformation mechanisms that 
are thought to operate only in very massive haloes. Fur- 
ther su pport for this was pr o vided bv Ivan den Bosch et al.l 
(|2008bl ) and iPasguah et al.l (|2008l ). who showed that, at 
fixed stellar mass, the average colours and concentrations, 
as well as the occurrence of star formation and AGN activ- 
ity, reveaJ_only_a^_\^y_w^^ dependence on halo mass (but 
see lWeinmann et al.ll200^. 

IWeinmann et all (|2006bl ) compared the fractions of 
red and blue galaxies in the SDSS group catalogue of 
IWeinmann et al.l ll2006ari with the s emi-a nalytic model of 
galaxy formation of ICroton et al.l (|2006l ). Although this 
model accurately fits the global statistics of the galaxy pop- 
ulation, the model predicts a red fracti on of satellites that 
is much higher than observed (see also iBaldrv et aLll2006l : 
ICoil et al.ll2003 ). 

In this paper we extend the study of IWeinmann et al.l 
(|2006bl ) using the much larger galaxy group catalogue of Y07 
and a larger suite of semi-analytic models. Another new as- 
pect of our work here is that we augment the SDSS-based 
galaxy properties in the group catalogue with information 
derived from the GALEX-SDSS matched sample of Salim et 



al. (2007; hereafter S07). The S07 analysis provides comple- 
mentary quantities such as stellar masses and star formation 
rates based on the two GALEX UV bands plus five-band 
SDSS photometry. The UV provides a much more sensitive 
probe of the recent star formation history of a galaxy than 
optical colours alone, which refiect the highly degenerate ef- 
fects of stellar populations, metaU icity, and dust extinction 
(|Yi et al.ll2005l : lKavirai et al.ll2007l ') . In addition, the UV can 
provide reliable measures of SFR for galaxies with weak or 
undetected emission lines. 

The goal of this paper is to investigate the dependence 
of the fraction of quenched galaxies on galaxy properties and 
DM halo mass in both the observational group catalogue and 
in several different semi-analytic galaxy formation models. 
Another new asp ect of our study with re spect to previous 
comparisons (e.g. IWeinmann et allbOOGbl ) is that we com- 
pare with the results from three independently developed 
semi-analytic codes, and for one of the codes, we examine 
three different variants with different physical ingredients. 
In this way, we investigate how these empirical results can 
constrain the input physics in these kinds of models. An out- 
line of the rest of our paper is as follows: in §2, we describe 
the empirical data sets and the group catalogues used in 
our study; in §3 we describe the theoretical models; in §4 we 
present the results of our comparisons between the empirical 
data and the models; we discuss our results and conclude in 
§5. 

2 THE DATA SAMPLE 

2.1 The SDSS Group Catalog 

We make use of the Galaxy Group Catalog of Yang et 
al. (2007; Y07). The catalogue was c onstructed by ap - 
plying the halo-based group finder of lYang et al] (120051 ) 
to the New Y ork University Value -Added Galaxy Catalog 
(NYU-VAGC. iBlanton et al.ll2005al l. which is based on the 
Sloan Digital Sky Survey (SDS S) Data Release 4 (DR4; 
lAdelman-McCarthv et all l2006h . From the Main Galaxy 
Sample, Y07 selected galaxies with extinction-corrected r- 
band apparent magnitude brighter than r = 17.77, within 
a redshift range 0.01 < z < 0.2, and with a redshift com- 
pleteness C > 0.7. They augmented this sample with 7091 
galaxies with 0.01 < z < 0.2 with redshifts from alternate 
sources. The resulting sample (Sample II of Y07) has a mean 
redshift of 2: ~ 0.1 and a total sky coverage of 4,514 deg^, 
and contains 369,447 galaxies. In this paper, we refer to this 
as "the SDSS sample". 

All absolute magnitudes are k and evolution (e) 
corrected to the z = 0.1 rest-frame, as described by 
iBlanton et"al] (|2003l l. Stellar masses are computed using the 
relation between rest-fram e optical colour and mass-to-light 
ratio of lBell et ah! (|2003l '). as specified in Y07, assuming a 
Kroupa IMF. Our sample is not volume limited, and we 
therefore attempt to correct for the resulting Malmquist 
bias by weighting e ach galaxy by a standard V^ax correc- 
tion (|Schmidtlll96i ). 

The group finder first identifies a p otential group cen - 
tre by the friends-of- friends algorithm (|Davis et al.1119851 ). 
and then computes the characteristic luminosity of each 
tentative group. The characteristic group luminosity is de- 
fined as the (incompleteness corrected; see Y07) combined 
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luminosity of all group members with "'^Mr ^ —19.5 + 
5 log h. The characteristic group stellar mass is similarly 
the incompleteness-corrected total stellar mass contributed 
by galaxies with °'^Mr ^ —19.5 + 5 log ft. Then, the veloc- 
ity dispersion and the virial radius of the dark matter halo 
(DMH) associated with each tentative group are estimated 
iteratively, assuming a constant mass-to-light ratio as an ini- 
tial guess. T hese two estimates can be used fo r determining 
the spherical iNavarroT Frenk. fc Wliit3 (|l997l ) (NFW)-type 
dark matter profile. Assuming a Gaussian probability distri- 
bution along the redshift direction and a projected spherical 
NFW dark matter profile for the perpendicular plane, group 
members are updated until there is no change in their mem- 
berships. 

Lastly, a DMH mass is assigned to each group, assum- 
ing a one-to-one correlation between the characteristic stel- 
lar mass for each group and t he halo mass derive d from a 
theoretical mass function (e.g. I Warren et al]|2005l ). We use 
the halo mass estimated from stellar mass, as luminosity- 
based halo mass estimates may be biased differently for 
groups with blue vs. red galaxies (Y07). Y07 find that they 
can assign group masses down to a lower limit of Mhaio < 
10^^ ®Mq//i, although the completeness of the group cata- 
logue begins to drop below unity for halo masses less than 
Mhaio ^ 10'^^ft~'^MQ. The central galaxy in each group is 
identified as the galaxy with the largest stellar mass. 

There are 204,813 groups in the resulting group cat- 
alogue. Note that we call all clusters "groups" regardless 
of their richness or density, including groups that contain 
only a single member. For a more detailed description of 
the group finder, and the results of extensive tests of its 
completen ess, contamination , and purity, we refer interested 
readers to lYang et all l|200iD . 



2.2 GALEX 

The GALEX data provide near-UV (NUV, effective wave- 
lengt h ~ 227lA) and far -UV (FUV, ~ 1528A) band informa- 
tion (|Martin et al.ll2005l '). Of 741 sq. deg. of GALEX unique 
imaging, 645 sq. deg. overlaps with the SDSS DR4 spec- 
troscopic area. We make use of the GALEX-SDSS matched 
sample constructed by Salim et al. (2007; S07). S07 first de- 
fine a SDSS parent sample in the GALEX overlap region 
of objects spectroscopically classified as galaxies, and hav- 
ing optical magnitude 14.5 < Vpstro < 17.77 and redshift 
0.005 < z < 0.22. This sample contains 49,346 galaxies. 
For each of the galaxies in the SDSS parent sample, S07 
searched within a 4 arcsec radius (corresponding to 7 kpc 
at z ~0.1) for a match in the GALEX source catalog. 

Using a large library of model SEDs based on the 
iBruzual fc CharlotI (|2003l ) population synthesis code, S07 
then performed Bayesian SED fitting to the 7-band photom- 
etry (FUV, NUV, u, g,r,i, z) to obtain estimates of dust ex- 
tinction, stellar mass, and star formation rate (SFR). These 
parameters were built from probability distribution func- 
tions, thus taking into account parameter degeneracies. The 
typical error in the specific SFR (obtained from the width 
of the SED fitting probability distribution) is 0.2 dex (star- 
forming galaxies) to 0.7 dex (passive galaxies). 

Of the full S07 sample, we use 32,787 galaxies that are 
associated with SDSS galaxies in our group catalog, and 
exclude a small number of objects with very poor fits (see 



discussion in §4.3 in S07). The redshift range for this final 
sample is 0.01 ^ z ^ 0.2 with mean redshift 0.104, and 
the redshift distribution is very similar to that of the parent 
SDSS sample. 



3 THEORETICAL MODELS 

We adopt a semi-analytic approach to model galaxy for- 
mation within the ACDM picture. We make use of a 
total of 5 sets of models based on different prescrip- 
tions for our comparison with the empirical data. Three 
of the models are constru cted using the latest ver- 
sion of the Somerville code feomor ville fc Primackl 19991 : 
ISomerville. Primack. fc Faber 2001; So merville et al.ll2008bl. 
here after S08). The others are the Millenn ium mod- 
els (jCroton et all |2006|; iDe Lucia et al.l 120061) and the 
MORGANA models (iMonaco. Fontanot. fc Taffonil l2007l : 
iFontanot et al.ll2006l . l2007l ). 

We describe the basic scheme in the Somerville code. 
We use an N-body simulation box to obtain the masses 
and positions of the "root" dark matter halos, and com- 
pu te the merger histor y for each halo using the method 
of ISomerville fc KolattI (|l999D . Within each dark matter 
halo, gas cools via atomic cooling (|White fc FrenkI Il99ll : 
ISomerville fc Primackl Il999l ) and forms a rotationally sup- 
ported disk. The r adial sizes of disks are com puted using the 
model described in ISomerville et al.l (12008^ 1 ■ which accounts 
for the initial Navarro-Frenk- White profiles of the DM halos 
and the "adiabatic contraction" due to the self-gravity of the 
infalling baryons. Cold gas is turned into st ars in the galactic 
disk following the Schmidt-Kennicutt law (|Kennicuttlll989l . 
Il998l ). and gas with s urface density lower than a critical 
threshold density (e.g. iMartin fc Kennicuttll200ll ) does not 
form stars. Massive stars explode as supernovae and reheat 
the cold gas. We trace chemical evolution using a simple "ef- 
fective yield" parameter. Each generation of stars produces 
a fixed "yield" of metals, which are deposited in the cold 
gas. This gas may then be ejected and mixed with the hot 
component by supernova or AGN driven winds. 

When a satellite is subsumed into a larger dark matter 
halo, it is assumed that it immediately loses its hot gas halo, 
and thus does not receive any new supply of cold gas. The 
orbital decay and eventual merging, due to dynamical fric- 
tion, of satellite galaxies within dark matter halos is tra cked 
using a modified version of the IChandrasekhail (| 19431 ) for- 
mula (|Bovlan-Kolchin. Ma. fc Quataert 2008h . Mass loss 
and tidal destruction are also accounted f or, using a sim- 
plified vers ion of the approach presen ted in iTavlor fc Babull 
j2004l ) and lZentner fc Bullockl (|2003l '). Galaxy mergers trig- 
ger bursts of star formation, the efficiency and timescales of 
which are modelled using results from hydrodynamic simu- 
lations of galaxy-galaxy mergers. 

The code also tracks the growth of black holes and the 
energy they produce. Every top level DM halo is seeded with 
a black hole of IOOM0. Mergers trigger the "bright mode" 
black hole growth that is associated with luminous quasars. 
Following every merger with mass ratio greater than 1:10, 
the black hole grows at its Eddington rate until it reaches 
a critical mass. The critical mass is that at which the en- 
ergy radiated by the black hole is sufficient to halt fur- 
ther accretion, i.e., the black hole growth is self-regulated 
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Table 1. Model descriptions 



Model 


box size (Mpc/h) 


Characteristics 


SOS AGN-FB (fiducial) 


120 


Bright mode AGN-driven winds + Radio mode AGN heating 


SOS no AGN-FB 


120 


Control model without AGN feedback 


808 HQ 


120 


cooling quenched according to the halo mass 


dL06 


120 


Radio mode AGN heating 


MORGANA 


144 


Bright mode AGN-driven winds + Radio mode AGN heating 



iHopkins et al.ll2007l '). Soon after reaching the critical mass, 
the black hole enters a "blowout" phase, resulting in a de- 
cline in the accretion rate. The associated radiation can also 
drive winds that remove cold gas from the galaxy (see SOS 
for details). 

The SOS models also incorporate "radio mode" feed- 
back associated with low-efficiency accretion. The accretion 
rat e is computed us i ng th e isothermal Bondi flow model 
of iNulsen fc FabianI (|2000l '). In the presence of a quasi- 
hydrostatic shock-heated gas halo, the energy from this ac- 
cretion is assumed to drive radio jets that can heat the gas 
and partially or completely offset the cooling flow. 

The resulting star formation and enrichment histo- 
rie s are convolved with th e stellar population models 
of lBruzual~fc" CharlotI (|2003l ') to compute magnitudes and 
colours. We have adopted a Chabrier IMF. The impact of 
dust extinction is modelle d using an analytic model, as in 
iDe Lucia fc BlaizotI (120071 ). 

The SOS fiducial model includes all of these mechanisms. 
We also consider a "no AGN feedback model", which does 
not include either the bright mode or radio mode AGN 
feedback mechanisms, and the "Halo Quenching" model, in 
which cooling is shut off when a halo grows more massive 
than ~ 10^^MQ/h~^ (see Table 1 for a summary of all mod- 
els). The halo quenching model is included as an illustration 
of a quenching mechanism that has a simple dependence 
on halo mass only. A similar model has been considered by 
ICattaneo et a l.' ('2006'). It i s bas ed on the ideas propose d 
bv lBirnboimfc Dckcl ()2003l ') and lDekel fc BirnboimI (|2006h . 
who suggest that whenever a halo grows above this critical 
quenching mass, the gas is shock-heated to near the virial 
temperature, and can be easily kept hot either by an AGN 
or by other pr ocesses such as heating by gas clumps or orbit- 
ing s atellites (|Dekel fc Birnboim|[2008l : iKhochfar fc Ostrikeij 
120081 '). 

We also consider two additional models from other 
groups in this study. The de Lucia et al. (2006, hereafter 
dL06) models contain similar ingredients to the SOS models, 
with the following differences. They are based on merger 
trees extracted from the Millennium N-body simulations 
(|Springel et al.ll2005l ). Unlike in the models of SOS, they 
do not include the effects of adiabatic contraction and an 
NFW halo profile in their estimates of galaxy sizes, which 
leads to a different evolution in the star formation rates 
and gas fractions in their models. They use somewhat dif- 
ferent (though similar in spirit) recipes for star formation 
and supernova feedback, and they use a different approach 
for modelling black hole growth (though, like SOS, they as- 
sume that "bright mode" black hole growth is triggered by 
mergers) . They do not include "bright mode" AGN feedback 
(AGN-driven winds) , and their treatment of the radio mode 



feedback is again similar in spirit but different in detail from 
SOS. Magnitudes and colours (including dust extinction) are 
computed in a similar manner to SOS, and use a Chabrier 
IMF. We obtained the dL06 catalogs from the public Mil- 
lennium database ( http: / /www. g- vo.org/Millennium) . 

We also consider the predictions of the semi- 



analy t ic model MORGANA ( Monaco. Fontanot. fc Taffonil 
I2OO7I : iFontanot et~al] |2006| . 120071 '). MORGANA follows a 
scheme similar to SOS and dL06, but it includes a differ- 
ent treatment for the thermal processes acting on baryonic 
gas. More details on the updated version we use in this paper 
are presented in Lo Faro et al. (2008, in prep). The model 
is based on merger trees obtained using the pinocchio 
method (|Monaco. Theuns. fc Taffonill2002 ) , similar to those 
predicted by N-body simulations. Gas coolin g and infall 
follow the prescription described and tested in IViola et al.l 
(|2008l ). while star formation and stellar feedback are then 
modeled as in 

Im 

onac 3 (|2004 ). When two DM halos merge, 
dynamical friction, tidal stripping and tidal shocks on the 
satellite galaxies are followed using the lTaffoni et al.l (|2003l ') 
formulation. Similarly to SOS, when a satellite DM halo 
merges into a larger one, all of its hot gas is shock heated 
according to the new halo potential and gets removed from 
the satellite (thus implying that the corresponding satellite 
galaxy does not receive any further col d gas supply). Disc 
sizes are computed using the model of IMo. Mao fc Whitg 
(|l998h . and bulge sizes are computed assuming that kinetic 
energy is conserved in mergers. 

A key ingredient in MORGANA is the inclusion of a 
self-consistent model for the accretion of gas onto su- 
permassive black holes and the resu lting AGN feedback 
mo des (following lUmernural 1200 ll and iGranato et al. I I2OO4 
see iFontanot et alT 20061 for more details) . This modeling 
assumes that the loss of angular momentum is the main 
regulator of black hole accretion. This is triggered by the 
presence of gas in the bulge component, and the almost com- 
plete loss of angular momentum of accreted gas is related 
to star formation activity. Following iGranato et al. 1 (|2004 ). 
star formation creates a reservoir of low angular momentum 
gas which is then accreted at a rate regulated by the vis- 
cous accretion time scale or by the Eddington limit. The 
nature of feedback from the AGN depends on the accretion 
rate in Eddington units; whenever this is higher than 0.01 
("bright mode") , the AGN can trigger a m assive galactic- 
scale wind (see [Monaco fc Fontanotl l2005l ) which leads to 
the complete removal of the ISM from the galaxy, while at 
lower accretion rates ("radio mode") the energy is ejected 
through jets that feed back on the hot halo gas with an ef- 
ficiency that scales with V^, where Vc is the halo circular 
velocity. As a consequence, BH accretion requires some star 
formation to be triggered, and feedback follows the onset 
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Figure 1. The distribution of ''■^r-band absolute magnitude vs. "'^(g — r) colour. Gray shading and green contours show the conditional 
probability PC'^A/rl^'^Cff — r)). The orange line shows the demarcation line for the "red" and "blue" galaxy populations used in this 
study. 



of cooling only after some time. The ejected energy heats 
the hot halo gas component and quenches the cooling flow. 
Galaxy SEDs, magnitudes and colours are obtained using 
the GRASIL spectro-photometric co de with radiative trans- 
fer for computing the effect of dust (jSilva et al.lll99d ). 

In Table 1, we present a brief description of each model. 



4 RESULTS 

In this section, we investigate the joint dependence of star 
formation quenching on stellar mass and halo mass, in or- 
der to try to contrain the physical mechanisms that are re- 
sponsible for quenching. We make use of two indicators of 
quenched star formation: red optical colours, and low spe- 



cific star formation rates {SSFR = SFR/Mga_i). Red optical 
colours are frequently used to isolate "quenched" galaxies; 
however, a red optical colour can arise from a degenerate 
combination of an old stellar population, a high metallicity, 
or strong dust extinction. SSFR based on UV-optical data 
are a more sensitive probe of recent star formation. 

In order to mimic the selection effects of the flux- 
limited observational sample, we first assign redshifts to 
all the model galaxies by placing an "observer" in a cor- 
ner of the simulation box. We apply the same fiux limit 
used in the observational catalogs to the models by selecting 
only galaxies with apparent r-band magnitude r < 17.77. 
We then apply the Vmax weighting factor, just as we do 
with the galaxies in the observational sample. Then, we ex- 
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Figure 2. Galaxy stellar mass vs. specific star formation rate. Gray shading and pink contours show the conditional probability 
P{Mg^l\SSFR). The purple line shows the demarcation line between "passive" and "active" galaxies used in this study. 



elude all halos that do not contain any galaxy brighter than 
"'^Mr ^ —19.5+5 log h, as these halos would not be included 
in the group catalog. 

We apply these selection criteria for all of the model- 
data comparisons shown in the main text, unless otherwise 
noted. We also present our main results for the models with- 
out these selection criteria in the Appendix. 



4.1 Global Distribution Functions: colour and 
SSFR 

In Fig. [T] we present the global colour-magnitude relations 
(CMRs) at z = 0.1 for the observations and theoretical mod- 
els, along with the dividing line between the observed red 
sequence and the blue cloud (sometimes called the "green 



valley"). In this figure, we show the results for the whole 
SDSS sample, regardless of inclusion in the group catalog, 
and similarly we have not applied the observational selec- 
tion criteria to the theoretical models. The magnitudes and 
col ours are shown in th e rest-frame z = 0.1 system defined 
bv lBlanton et al.l (|2003h . For the dL06 models, we used the 
observed frame magnitudes at z = 0.1, converted to absolute 
magnitudes. For the MORGANA galaxies, we compute the ab- 
solute z = 0.1 magnitudes from the corresponding synthetic 
spectra. For the SOS models, when we computed standard 
z = 0.0 frame colours and magnitudes, we produced a good 
match to the observed colour-magnitude relation expressed 
in the z = frame. However, when we computed the z = 0.1 
system colours as described in lBlanton fc RoweisI (|2007l ) . we 
found it necessary to apply a shift of 0.05 magnitudes to the 
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"'^{g — r) colour to match the location of the observed z = 0.1 
system red sequence. This is indicative of small differences in 
the shape of the SED's i n the semi-ana l ytic r nodels from the 
synthetic SED's used bv lBlanton et al.1 1\200^ for computing 
the k-corrections. It should also be noted that details in the 
population synthesis prescriptions, such as chemical enrich- 
ment and dust extinction, may also cause noticeable differ- 
ences in colours (see Appendix for the effect of dust). There- 
fore, reproducing the CMR quantitatively in semi-analytic 
models is quite challenging. 

The observations clearly show the familiar red sequence 
and blue cloud, and these features are reproduced rea- 
sonably well in all of the theoretical models, except the 
SOS no-A GN-feedback model. As has been pointed out be- 
fore (e.g., ICroton et al] l2006l : ICattaneo et al.1 120061 ) ■ semi- 
analytic models without some kind of suppression of cooling 
in massive halos, e.g. by AGN feedback, predict that massive 
galaxies are still accreting plenty of cold gas at the present 
day, and therefore are star-forming and blue, in conflict with 
observations. We do see subtle differences in the structure 
of the CM distribution for the different models, for exam- 
ple, the S08-halo-quenching model produces very few low- 
luminosity red galaxies, and the dL06 model produces a very 
strong bimodality. Also, we notice that the red sequence in 
the halo-quenching model is slightly bluer compared with 
the SOS fiducial model. This is because the central galaxies 
in the halo-quenching model stop forming stars once their 
halo becomes more massive as 10^^ Mq. In general, this 
happens at an earlier epoch than the "radio mode" feedback 
is able to quench star formation. Because chemical evolution 
is also halted when star formation is quenched, the massive 
galaxies in the halo-quenching model do not become as en- 
riched as do the galaxies in the fiducial model (see SOS). 
We could have adjusted for this by increasing the chemical 
yield, but we chose to leave all the free parameters fixed in 
both models to allow a direct comparison. 

The behavior of MORGANA is somewhat different from 
that of the other models. A clear red sequence is present 
at intermediate and faint magnitudes, but it fades away at 
bright magnitudes, in contrast to the observations. Lumi- 
nous galaxies tend to be blue. We interpret this result as 
an inefficient quenching of star formation through "radio- 
mode" feedback in morgana. In fact, as explained in §. |3l 
in this model the AGN heating switches on only after some 
cooled gas has already started forming stars in the host 
galaxy. Obviously the residual activity is stronger for longer 
time delays between the onset of cooling flows and the accre- 
tion onto the central black hole. On the other hand, the blue 
cloud seems depleted. As shown in Fontanot et al. (200S, in 
prep), the depletion of the blue cloud is connected in part 
to the "satellite overquenching problem", discussed further 
in §5, caused by the too-efficient strangulation of satellite 
galaxies. However, we shall see later that MORGANA also pro- 
duces too few low mass blue central galaxies. This is likely 
due to strong supernova feedback. 

Given the strongly bimodal colour distribution exhib- 
ited by both the empirical data and most of the mod- 
els, it is natural to define a dividing line between the two 
populations, and to investigate the fraction of galaxies in 
the red population, /red, as representing quenched objects. 
We adopt the same demarcation line between red and blue 



galaxies as Weinmann et al. (2006a, hereafter W06a): 

°'^{9~r) = 0.7 - 0.032('''^A/r-51ogio/i-f 16.5) (1) 

as shown in Fig. [T] 

We can also define the specific star formation rate 
(SSFR), as the present star formation rate divided by the 
stellar mass of the galaxy (= SFR/Mg^i), and plot a simi- 
lar diagram in terms of SSFR and mass. In Fig. (2] we show 
the conditional probability distribution for SSFR as a func- 
tion of stellar mass, as derived from the GALEX+SDSS ob- 
servations, and for the models. The star formation rate is 
averaged over the past 100 Myr. Note that the number of 
galaxies in the sample is roughly S times smaller than the 
sample shown in the colour-magnitude distribution plot, be- 
cause the field coverage of GALEX survey used in the Salim 
et al. study is not as wide as that of the SDSS. In a similar 
manner, we define a cut in SSFR to separate "active" star 
forming galaxies from "passive" ones. We define a galaxy as 
"passive" if the following condition is met: 

logio SSFR < -9 - 0.2 log[Mg^i /h^^Ms] (2) 

Note here that our criterion is based on galaxy stel- 
lar mass in order to make a more direct comparison with 
the theoretical models. It is also worth noting that our cri- 
terion roughly corresponds to SFR ~ lM0yr~'^ at Afgai'^ 
IO^/i-^Mq and SFR ~ ClMgyr"^ at Mgai- IO^-'^/i-^Mq. 
Below O.lMQyr"^, SFR obtained from multiband photom- 
etry may not be robust due to the degeneracy of burst 
time and the mas s fraction of a young population (e.g. 
iKavirai et allbOOTt ). 

Our demarcation is comparable to that of W06a when 
their criterion, which is based on the "'^r-band magnitude, is 
lowered by 0.6 dex. However, SFR estimates from emission 
lines, which were used by W06a, trace only very recent star 
formation, whereas the UV used here traces star formation 
over a longer timescale ('^ 1 Gyr). Therefore our "passive" 
sample is not directly comparable to that of W06a, though 
the results are qualitatively similar. The amount of UV flux 
in the passive galaxies is very small and can still be con- 
sistent with the amount of UV flux that can be produced 
by old stars s uch as low-mass horizontal-branch stars (see 
lYi et al]|2005l for details). In this regard it is justifiable to 
call them "passive". 

In Fig. (2] we see a well-defined "star-forming sequence" 
that produces the blue sequence in the traditional CM di- 
agram, but the "quenched" population that produces the 
tight red sequence is quite spread out in SSFR. This is sim- 
ply a reflection of the relative insensitivity of optical colours 
to small amounts of recent star formation. The precise values 
of SSFR at low levels of star formation cannot be estimated 
very accurately from the data, and therefore one should not 
take too seriously the position of the quenched galaxies in 
the SSFR vs. Mgai plot. Again, the models with AGN feed- 
back (or halo quenching) qualitatively reproduce the em- 
pirical distribution reasonably well, but show some inter- 
esting differences. The dL06 model has a strongly-bimodal 
distribution of SSFR, with most galaxies living either on 
the star-forming sequence or being completely quenched. In 
contrast, the SOS model has a larger population of "semi- 
quenched" galaxies at intermediate values of SSFR. The 
MORGANA model produces an even broader distribution of 
SSFR, with many quenched low-mass galaxies. 
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Figure 3. The fraction of red galaxies (fred; left) and passive galaxies (fpassive; right) as a function of galaxy stellar mass. The empirical 
data are shown by a black solid line, and coloured lines show the results of the theoretical models, as indicated on the plot. We show 
the dependence for all galaxies (top), and for central (middle), and satellite (bottom) galaxies seperately. Each point contains at least 
10 galaxies. Most of the SAMs reproduce the trend for central galaxies reasonably well, but predict a much larger fraction of small-mass 
red/passive satellite galaxies than are observed. The error bars indicate Poisson errors. 



4.2 Dependence of Star Formation Quenching on 
Stellar Mass 

Fig. |3] shows the dependence of /red and /passive on galaxy 
stellar mass. We show /rod and /passive vs. stellar mass for all 
galaxies (upper panels), and for central and satellite galax- 
ies separately (middle and lower panels, respectively). Al- 
though we do not specifically use the information from the 
group catalog in the uppermost panel, in all panels we use 
only the SDSS galaxies that are included in the group cat- 
alog, and apply the group catalog-like selection criteria to 
the theoretical models, as described at the beginning of this 
section. 

In the empirical data, /rod shows a strong dependence 
on galaxy mass for both centrals and satellites, in the well- 
known sense that low-mass galaxies are largely blue, whereas 
massive galaxies are more likely to be red. It appears that 
/red is a steeper function of stellar mass for central galax- 
ies than for satellites. In general, the red galaxy fraction 
is higher for satellit e s for a fixed galaxy mass. Similarly, 
Ivan den Bosch et al.l (|2008al lbh found that satellite galaxies 
have redder mean colours than centrals at a fixed stellar 
mass. The trends appear qualitatively very similar when we 
consider /passive as a function of stellar mass. 

The SOS fiducial model, SOS halo-quenching model, and 
dL06 model all reproduce the trends in /rod and /passive with 
stellar mass quite well for central galaxies. Note the similar- 



ity of the predictions of the SOS fiducial and halo-quenching 
model. The SOS no AGN feedback model, as we have already 
seen, predicts the presence of too many massive blue central 
galaxies. MORGANA reproduces the sense of the trend of the 
red fraction with stellar mass, but slightly overproduces red 
galaxies at low stellar masses, and significantly overproduces 
blue galaxies at high stellar masses. 

All of the models badly overproduce the number of low 
mass red satellites, and predict too weak a trend of /red with 
stellar mass for these objects. The dL06 models predict a 
nearly fiat run of /rod with stellar mass for satellites, with 
values that are much too high (close to unity) compared 
with the empirical values. At the highest satellite masses, 
there is too high a fraction of blue galaxies in the MORGANA 
model and the SOS no-AGN-feedback model. We can see by 
comparing the SOS no-AGN-feedback model with the fiducial 
model that AGN feedback does not affect galaxy colours 
below stellar masses of ~ lO^°/i~^M0; therefore the excess of 
low-mass red satellites is not connected with AGN feedback. 

The conclusions we would draw from the comparison 
with /passive are qualitatively similar, though quantitatively 
somewhat different. For example, the dL06 model shows a 
better quantitative match to the /passive data than to /red for 
central galaxies. The SOS fiducial and halo-quenching models 
produced almost indistinguishable results for /rod(Afgai), but 
significantly different results for /passive (Mgai). This is due 
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Figure 4. The red fraction f^od (left) and passive fraction fpassivc (right) as a function of halo mass, for different stellar mass bins, as 
indicated by different colours and symbol sizes (see plot legend). In each plot, results are shown for the observed SDSS or GALEX+SDSS 
group catalog (top left) and the five models as indicated on the plot. We present the results for central {top set of panels) and satellite 
{bottom set of panels) galaxies seperately. Each point contains at least 5 galaxies. The SAMs reproduce the main trends reasonably well 
for central galaxies, but satellite galaxies do not show the correct trend of /^cd with stellar mass. 



to the age-metallicity degeneracy — as we discussed in 
massive galazxies in the halo-quenching model are more metal 
poor than in the fiducial model. Therefore, although they are 
older (as seen in the /passive diagram), their optical colours 
are similar. Similarly, the large difference between /red and 
/passive for the SOS no AGN-feedback model is due to dust 
extinction: in this model, massive galaxies are actively star 
forming, and therefore extremely dusty (see Appendix for 
more discussion of the effects of dust extinction) . 

The somewhat lower fraction of red/passive satellite 
galaxies in the SOS models is due to the inclusion of tidal 
destruction, which is not included in the dL06 or IVIORGANA 
models. In the SOS model, satellites that orbit within their 



host halo for a long time can eventually become tidally de- 
stroyed, and their stars are added to a "Diffuse Stellar Halo" 
(see SOS). Naturally, in the absence of tidal destruction, 
these satellites exhaust all of their gas and become very red. 
However, we see here that although the inclusion of tidal de- 
struction helps reduce the excess of low-mass red satellites 
in the models, a significant discrepancy still remains. 

One might then wonder whether increasing the effi- 
ciency of tidal destruction could completely solve the satel- 
lite 'over-quenching' problem that we see here. We do not 
believe that this is a viable solution, for several reasons. 
First, the tidal disruption model used by SOS was tuned 
to match the sub-structure mass function for very high- 
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Figure 5. The red fraction fjcj (left) and passive fraction fpassivc (right) as a function of stellar mass, for different halo mass bins, as 
indicated by different colours and symbol sizes (see plot legend). The rest of the plot details are as in Fig. |4] 



resolution N-body simulations. SOS also showed that their 
models correctly reproduce the total number of low-mass 
galaxies (i.e. the faint end slope of the stellar mass func- 
tion), although the model overproduces low-mass bulge- 
dominated (red/passive) galaxies and underproduces low- 
mass disc-dominated (blue/active) galaxies. Significantly in- 
creasing the efficiency of tidal destruction would be in con- 
flict with the N-body results and would also produce an 
overall deficit of low-mass galaxies. Put another way, tidal 
destruction can remove low-mass red satellites but cannot 
increase the number of low-mass blue satellites. 



4.3 Dependence on Stellar Mass and Halo Mass 

In this section we explore the dependence of /red and /passive 
on DM halo mass and stellar mass. In Fig. |4]we show /red 
and /passive as a function of halo mass, for different bins in 
stellar mass. We show the results for /red and /passive for 
central galaxies (top row), and for satellite galaxies (bottom 
row). In Fig. (5] we show a similar plot, but this time with 
the galaxy stellar mass plotted on the x-axis, and different 
bins in halo mass shown by the different colors. 

In the empirical data, we see that for massive Mgai > 
lO^^/i^^M0 central galaxies, there is no significant depen- 
dence of /red on halo mass (environment) for fixed stellar 
mass. For intermediate mass galaxies (lO^^/i"^ < Mgai < 
lO^^/i~'^M0), the dependence on halo mass appears to be 
stronger, but there is a limited region of overlap in galaxies 
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with different stellar masses that occupy halos of the same 
mass. This is because there is a fairly tight correlation be- 
tween halo mass and stellar mass. Similar results are again 
obtained for /passive- 

The SOS fiducial and dL06 models, both of which 
include AGN feedback, do reasonably well at repro- 
ducing the overall trends for central galaxies, as does 
the SOS halo-quenching model. We see a hint of a de- 
pendence on stellar mass in intermediate mass halos 
(lO"'^/i~^M0 < Mhaio ^ lO^^-^ft-^Mo) in the SOS models, 
while in the dL06 model, the dependence seems to be al- 
most solely on halo mass. However, interestingly, we see al- 
most the same stellar mass dependence in the SOS fiducial 
and halo-quenching models, while we know that the quench- 
ing mechanism is a pure function of halo mass in the halo- 
quenching model. The SOS no- AGN feedback model shows 
the correct trend with stellar mass at fixed halo mass (more 
massive galaxies have higher /red), but the opposite trend 
with halo mass (more massive halos have lower /red). Inter- 
estingly, MORGANA shows similar trends to the SOS no-AGN- 
feedback model: this implies that the current implementa- 
tion of AGN feedback in this model is insufficient to fully 
cure the "star formation quenching" problem. 

Considering /passive, we see somewhat different behav- 
ior. The SOS no AGN-feedback model produces no quenched 
galaxies at all in terms of /passive- The SOS and dL06 mod- 
els appear surprisingly similar, both showing almost a pure 
halo mass dependence (no significant dependence on stellar 
mass). On the other hand, MORGANA predicts a significant 
dependence of /passive on stellar mass, while showing overall 

low values of /passive. 

Observed satellite galaxies show an /red dependence 
that is very similar to that of centrals for massive galax- 
ies (Mgai 10^°'^/i~^Mq). For intermediate and low mass 
satellites, /red shows some dependence on both galaxy stellar 
mass and halo mass, but does not show the very sharp drop 
over intermediate halo masses (11.5 Mhaio 12.5/i~^M0) 
seen in the central population. In the semi-analytic models, 
we see that /red for the satellites does not have a strong 
enough dependence on Mgai at fixed halo mass. 

The stellar mass dependence of /passive for satellites in 
the empirical sample is not as clear as it was in terms of 
/red. This may be in part due to the smaller size of the 
GALEX-SDSS matched sample used to obtain /passive. The 
SOS models now all show a clear inverted trend: /passive is 
higher for lower mass galaxies (the opposite of the empirical 
trend). For the dL06 model, nearly all satellites are passive 
regardless of their stellar mass or halo mass. In MORGANA, 
satellite properties are a weak function of stellar mass, and 
the values of /passive are overall too high. 

Based on Fig. |4] alone, we might be tempted to con- 
clude that quenching is primarily a function of halo mass. 
However, Fig. [5] shows that quenching could equally well 
be considered to be primarily a function of stellar mass. We 
conclude that, especially for central galaxies, the degeneracy 
between stellar mass and halo mass is too strong to reach a 
firm conclusion on this point. 

In Fig. |6l we present the /rod and /passive distributions 
in the (Mhaio, Mgai) plane. We pixelise the (Mhaio, Mgai) 
plane, compute /red and /passive in each pixel, and indi- 
cate its value by the colour of the pixel. For example, a 
red colour indicates that galaxies are mostly red or pas- 



sive within the pixel, while a blue colour indicates that 
most of the galaxies are blue/active. These diagrams re- 
veal a number of interesting features. Considering the dia- 
gram for central galaxies in the empirical sample, we see that 
above a critical halo mass (Mhaio ^ lO^^/i^^Mo), nearly all 
galaxies are red and passive. For intermediate halo masses 
10^^ ^ Mhaio ~ lO^^/i~^M0, the structures show a complex 
dependence on both halo mass and stellar mass. In par- 
ticular we note that above a "critical" stellar mass of 2 — 
3 10^''/i~^Mq, the majority of galaxies are red and passive, 
regardless of their halo mass (though such massive galax- 
ies are not found in halos less massive than 10^^/i"^Mq). 
Comparing with the models that showed good qualitative 
behavior in terms of the binned quantities (SOS fiducial, SOS 
halo-quenching, and dL06), we can see that the distribution 
of the patterns in (Afhaio, Mgai) space are quite different from 
the observations — in general, the structures show stronger 
vertical divisions, indicating a stronger dependence on halo 
mass than on galaxy mass. It is interesting to note that 
these three models look much more similar to one another 
than any of the models does to the empirical data. It is also 
interesting that in terms of /red, the SOS fiducial and SOS 
halo-quenching models look very similar to one another, but 
they look extremely different in the /passive diagram. This 
again illustrates that optical colours are not an ideal probe 
of star formation quenching. 

The observed satellites show an interesting striation, 
which is nearly horizontal at the highest and lowest stel- 
lar masses, but somewhat diagonal for intermediate masses. 
It appears that the majority of satellite galaxies are 
red/passive if they are more massive than a few times 
lO^''/i~^M0 (just as for the central population), and are 
predominantly blue/active if they are less massive than 
lO^''/i~^M0, regardless of their halo mass. For intermedi- 
ate halo masses, it seems that unlike for central galaxies, 
the critical mass that marks the transition between mostly 
blue and mostly red galaxies is a function of halo mass, and 
is lower for higher halo masses. This suggests that the star 
formation activity in the most massive satellites is regulated 
by the same processes that shape centrals, while lower mass 
satellites are infiuenced by environmental processes such as 
tidal forces or ram pressure stripping. 

All the models fail quite miserably to reproduce the 
satellite properties. In addition to simply predicting too high 
a fraction of red satellites, none of the models shows the 
diagonal pattern of contours in the /red or /passive diagrams. 
The few pixels with high blue fractions in the models lie at 
high stellar mass for their halo mass, which is the opposite 
from what is seen in the empirical data. These are simply 
galaxies that were forming stars as centrals and have become 
satellites very recently. 

Readers are referred to the Appendix for a discussion 
of the impact on our results of the modelling of dust extinc- 
tion, our imposed selection criteria, and possible biases in 
the procedure for assigning halo masses in the SDSS group 
catalog. 

4.4 Connection virith Morphology, Black Hole 
Formation and AGN Feedback 

These results naturally beg the question: which physical 
process(es) are responsible for imprinting this dependence 
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Figure 6. The fraction of "red" galaxies (/i-edi left) and passive galaxies (/passive; right) for central {top set of plots) and satellite [bottom 
set of plots) galaxies in the (Mjiaio, -'^gal) plane. The fraction of red/passive galaxies in a given pixel in (Mhaioi -'^-'^gal) is indicated by 
the colour, where red colours indicate a higher red/passive fraction, as shown in the scale. To guide the eye, we draw a solid (dotted) 
line showing the approximate upper envelope of the central (satellite) galaxy mass distribution for the observational group catalog, and 
repeat this same line on every panel. Central and satellite galaxies in the observational group catalogs show noticably different joint 
dependencies on stellar mass and halo mass. The models with AGN feedback qualitatively reproduce the trends for central galaxies, but 
do not reproduce the empirical trends for satellites. 



of star formation quenching on galaxy and halo mass? Al- 
though this is a complex question that we will not be able to 
fully address in this paper, we attempt to at least identify 
some promising hypotheses that can be pursued futher in 
the future. 

As already disc ussed by many authors (e.g. 
iKauffmann et all hooj ). we are suspicous that the corre- 
lation of star formation quenching with stellar mass may 
in fact be linked to the tendency of the increasing bulge 
fraction also with stellar mass and perhaps halo mass. 
Although we have only rough morphological information for 



SDSS, we make use of a standard cut in concentration index 
(C = rgo/rso) to coarsely divide our sample into early and 
late type galaxies (early type s have C > 2.6; IStrateva et al.l 
I2OOII : IShimasaku et alTl200ll '). Such classification using the 
concentration i ndex is subject to co ntamination at roughly 
the 20% level (jStrateva et al.ll200lh . We then plot /early in 
the (Mhaio,Mgai) plane as befor^E, and show the results in 



^ Because of the difficulty of mapping the morphological infor- 
mation available in the models to the observable concentration 
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Figure 7. The fraction of early- type galaxies (fcarly) for centrals 
{left) and satellites {right) in the SDSS group catalog. Early type 
galajcies are defined as having a concentration index (C) greater 
than 2.6. 



Fig. [71 We see a strikingly similar pattern to the one seen 
when we plotted /rod and /passive in this way in Fig. [B] 
indicating a very strong correlation between red, passive, 
and early type galaxies. 

We test the connection of black hole mass with star 
formation quenching by plotting the ratio of black hole mass 
to the total stellar mass of the galaxy ((MBjf/Mgai)) in the 
same (Mhaio,Mgai) plane, which we do in Fig. [8] (for central 
galaxies only). Note that all three semi-analytic model codes 
reproduce the empirical relation between bulge mass and 
black hole mass reasonably well within the observational 
errors. In the SOS, dL06, and MORGANA models, we see that 
the galaxies that have the highest stellar mass to halo mass 
ratios have the lowest black hole mass to stellar mass ratios. 
This is a reflection of the spread in halo merger histories 
at fixed halo mass. In halos in which a massive black hole 
is formed relatively early, the cooling flow is also shut off 
at an earlier time, halting further galaxy growth except by 
mergers. Conversely, halos with relatively small black holes 
for their mass will be able to continue to cool, and the central 
galaxies will continue to form stars and remain blue. This 
is seen to be the case in the "blue ridge" in Fig. [S] (see 
panels h, k, and 1). Note that neither the strong trend in 
{Mbh / Mgai) nor the "blue ridge" is seen in the SOS "halo 
quenching" model, in which quenching is regulated purely 
by the halo mass. 

In all three models with AGN feedback, there is also 
a trend with halo mass in the sense that larger mass ha- 
los have larger black hole-to-stellar mass ratios. This is due 
to the dual modes of black hole growth in the models. In 
the "bright mode", the growth of the stellar bulge and the 
black hole are linked. Most of the black hole growth in the 
models occurs via this bright mode of black hole feeding. 
However, at late times, large mass halos can develop a hot 
hydrostatic halo which is assumed to fuel the "radio mode" 
of black hole growth. In these halos, the black hole can grow 
without any associated star formation, leading to an increase 
in (Ms£f/A/gai). This is supported by the fact that we do 
not see such a trend in the halo quenching model, which 
only contains black hole growth via the bright mode. 

We attempt to investigate whether this trend exists 
in the empirical data, using the SDSS velocity disper- 



index, we do not show the model predictions. The qualitative 
trends in the models are similar. 
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Figure 8. Average black hole mass to galaxy stellar mass ratio 
{MsH/^gal) for central galaxies, shown by the colour coding, as 
a function of Mjiaio a-^d Mgai, for four of the theoretical models. 
For the SOS and dL06 models, note the similarity of the structures 
seen here to those seen in the plot of /passive in Fig- E] This sug- 
gests that for intermediate halo masses, /passive is closely related 
to {Mgjj /Mgai) in the models. We also show {Ma/Mgai) for the 
SDSS group catalogs, where Ma, based on the SDSS measured 
velocity dispersion and the observed Mbh — f relation, is used as 
a proxy for the average black hole mass. Interestingly, the strong 
dependence of (A/o-/A/gni) on halo mass is not visible. 



sion (T as a proxy for black hole mass. We compute Mo- 
fr om the Mbh — Oe rel ation using the empirical relation 
of iGebhardt et all (j2000l l. We assign a scatter to the black 
hole mass at a given a by choosing a uniform rand o m de - 
viate over the Icr range quoted by IGebhardt et al.l (|2000f ) 
(we obtain indistinguishable results when a Gaussian ran- 
dom deviate is used). Interestingly, we see no evidence of a 
trend in {Mbh /M^ai) with halo mass. 

However, this result is not strongly conclusive, because 
it is not known whether the relationship between Mbh and 
(Je depends on halo mass. However, these results are sugges- 
tive that there may be too much black hole growth via the 
"radio mode" in the models, perhaps indicating that other 
heating processes not associated with black hole growth may 
also be important in quenching cooling flows. 

We caution, as well, that although it is tempting to try 
to interpret the "fine features" in these diagrams, the as- 
signment of halo masses in the current SDSS group catalogs 
is not very precise (see the Appendix), and therefore only 
broad statistical trends should be taken seriously. It is possi- 
ble that in the future, if more precise estimates of individual 
halo masses become available for large samples, for exam- 
ple from gravitational lensing or X-rays, we may be able to 
investigate these predicted trends in more detail. 
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5 DISCUSSION AND CONCLUSIONS 

In this paper, we set out to investigate the significance of in- 
ternal galaxy properties vs. environment in shaping the star 
formation history of galaxies, and to attempt to understand 
some of the physical processes that might be at work. We 
made use of the observational Group C atalog constructe d 
from SDSS DR4 and the NYU-VAGC d^n^eTaO [20o3) • 
The Group Catalog provides an estimate of the halo mass 
for each group, which can be directly compared with semi- 
analytic models. These halo mass estimates should provide 
a more unbiased probe of global environment than measures 
based on localgalaxy density (e.g. distance to the n-th near- 
est neighbor) Q 

We have also used a sub-sample of S DSS with GALEX 
coverage to estimate star formation rates (|Salim et al.ll2007l ). 
which should be a more direct probe of the physics of 
star formation quenching than optical colours. Besides, we 
make use of semi-analytic models from several independent 
groups (SOS, dL06, MORGANA), and containing different sets 
of recipes representing physical processes. 

We first investigated the global distribution of colour 
vs. magnitude and SSFR vs. stellar mass in the five models 
we considered, compared with the empirical data. Although 
the different models showed some differences in the details 
of their predictions for these quantities, the models with 
some kind of quenching (either due to AGN feedback or 
according to a critical halo mass) have similar qualitative 
features and on the whole are a reasonable match to the 
empirical data — not surprisingly, as these observational 
quantities have been a target for theoretical models for some 
time. We therefore found that it was sensible to identify 
active or quenched galaxies according either to a colour- 
magnitude criterion (the usual "green valley"), or a similar 
SSFR-Mgai criterion. 

Next we investigated the stellar mass dependence of 
the quenched fraction, based on optical colours (/red), 
or on SSFR determined from SDSS+GALEX photometry 
(/passive). Here we largely confirmed and reproduced results 
shown previously by other authors (though our results in 
terms of SSFR from GALEX are new) , namely that the frac- 
tion of red/passive galaxies increases with stellar mass for 
both central and satellite galaxies, and that the models with 
AGN feedback or halo mass-based quenching reproduce this 
trend reasonably well for central galaxies, but fail badly for 
satellites. All of the models produce too high a fraction of red 
satellites and too fiat a dependence of quenching on stellar 
mass, which we term the satellite over-quenching problem. 

We then investigated the joint dependence of quenched 
fraction (/red and /passive) on galaxy mass and halo mass. 
First we investigated /red and /passive as a function of halo 
mass, in different stellar mass bins. A difficulty with this 
approach was that, especially for central galaxies, there is 
quite a limited range of stellar masses in halos of a given 
mass. Our analysis showed that, for central galaxies in the 
observational group catalogs, the fraction of quenched galax- 
ies shows a strong dependence on halo mass at fixed stellar 

^ But note that some groups have started measuring galaxy num- 
ber densities using search ellipsoids la rge enough to be repre- 
sentative of mas sive cluster halos (e.e... ISchawinski et al.ll2007al : 
lYoon et al.ll2008l) 



mass, but also shows a strong dependence on stellar mass at 
fixed halo mass. We were not able to determine which quan- 
tity is the primary driver of quenching. The SOS fiducial 
and dL06 models reproduced these trends fairly well. For 
observed satellite galaxies, there was a much stronger de- 
pendence on stellar mass visible in the /red diagram (based 
on optical colours), and a weaker trend in the /passive dia- 
gram (based on UV-derived SSFR). Once again, all models 
failed to reproduce the satellite properties, and even showed 
an inverted trend in /passive with respect to the empirical 
data. 

We also found it interesting to look at the pattern of 
/red and /passive iu terms of the two dimensional (Afhaio- 
Mgai) plane. This analysis revealed that the contours of /red 
and /passive for central galaxies can be interpreted either 
as a horizontal run or as a vertical run, again due to the 
halo mass - central galaxy mass degeneracy. The models 
showed quite a different pattern in this space, and tended 
to show stronger vertical boundaries, indicating a stronger 
dependence on halo mass. For both the empirical data and 
the models, these diagrams demonstrate the complexity of 
the interplay between halo mass and stellar mass, and are a 
promising tool for posing stringent tests on physical recipes 
in galaxy formation models. However, we caution that the 
estimates of halo mass in the SDSS group catalogs are sta- 
tistical in nature, and this may introduce distortions into 
these distributions (see Appendix). 

We attempted to probe the physical origin of these re- 
sults by exploring additional correlations, such as the frac- 
tion of morphologically early type (spheroid dominated) 
galaxies in the observed sample, /early We found a strik- 
ingly similar pattern for /early in the (Mhaio,Afgai) plane as 
we did for /red and /passive , suggesting that these quantities 
are tightly linked in some way. One natural possibility is that 
the bulge mass is correlated with the mass of a supermas- 
sive black hole, and that the black hole mass in turn controls 
the quenching of star formation. Intriguingly, we found that 
in the models with black hole-regulated AGN heating (SOS 
fiducial and dL06), the galaxies that were most likely to be 
blue, actively star forming, and disk dominated, were ex- 
pected to be those with the smallest black hole for their 
mass. In addition, we saw a trend with halo mass, in the 
sense that central galaxies in larger mass halos are able to 
grow black holes more efficiently (the ratio of black hole mass 
to stellar mass is larger) . These trends were much weaker or 
absent in models in which the black hole is not involved in 
regulating cooling (such as the "halo quenching" model). We 
do not have direct estimates of black hole masses for large 
samples of galaxies, but we used the measured SDSS veloc- 
ity dispersion Oe and the observed Mbh — o"e relation to 
obtain estimates of a black hole mass proxy. Mo-. We do not 
see a strong trend in Ma/M^^i for the SDSS group catalog, 
indicating either that this method for estimating the empir- 
ical black hole masses is too crude, or that the dependence 
of black hole mass on halo mass in the models is too strong. 

Although the model predictions for the distribution of 
/red and /passive iu the (Mhaio, Afgai) plauc do not match the 
observational results in detail, we conclude that the observa- 
tional data are consistent with the basic qualitative picture 
presented by the models in which cooling is regulated by 
AGN feedback at least for central galaxies. In these mod- 
els, the suppression of cooling and the quenching of star 
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formation depend on two factors: the presence of a quasi- 
hydrostatic hot halo (strongly correlated with halo mass), 
and the mass of the supermassive black hole (strongly cor- 
related with galaxy mass). This picture is also str ongly sup- 
ported by the recent work of lPasguali et al] (|2008l 'l. which di- 
rectly explored the dependence of "radio mode" and "bright 
mode" AGN activity on halo mass and stellar mass in these 
same SDSS group catalogs. The empirical data do not seem 
to support models in which the process t hat suppresses cool- 
ing is solely a function of halo m ass fe.e:.. lDekel fc BirnboimI 
l2008l : lKhochfar fc OstrikeilbOOSl 'l. although it will be impor- 
tant to explore the explicit predictions of alternate heating 
mechanisms (such as heating by clumps or infalling satel- 
lites) in detail. 

The MORGANA model suffers from the largest disagree- 
ment with the observations. The treatment of the triggering 
of "radio mode" accretion in MORGANA is significantly dif- 
ferent than in the other two models, requiring that star for- 
mation is inevitably associated with the triggering of radio 
mode accretion. Although the radio mode feedback mecha- 
nism adopted in MORGANA is largely able to solve the "over- 
cooling problem" in terms of reproducing the galaxy stel- 
lar mass or luminosity function, this star formation makes 
many massive galaxies too blue. Our analysis places tight 
constraints on the links between star formation and AGN 
activity at late times, and highlights our current lack of un- 
derstanding of the details of the processes that regulate both 
kinds of activity. 

For satellite galaxies, the empirical diagrams suggested 
that quenching depends on both stellar mass and halo mass, 
such that the "critical stellar mass" that divides active from 
passive (blue from red) galaxies decreases with increasing 
halo mass. None of the models was successful in predicting 
this trend. Including tidal destruction of satellites, as was 
done in the SOS models, improves the agreement with the 
data because satellites that have been orbiting for a long 
time within the host halo (which tend to be red and pas- 
sive) are destroyed. However, it seems that tidal destruction 
cannot provide a full solution to the problem. 

The problems with over-quenching of satellites in 
semi-analytic models h ave been demonstrated before (e.g. 
IWeinmann et aL 1 11' .Ij i and are likely due to the assump- 
tion applied in nearly all SAMs that the hot gas halo, 
which is the source of new cooling gas, is stripped off 
instantly when a galaxy becomes a satellite in a larger 
halo (sometimes called "strangulation"). Therefore satel- 
lites fairly quickly consume their remaini ng cold gas reser- 
voirs and become red and passive (e.g. ICrowl fc KennevI 
I2OO6I ). However, recent hydrodynamic simulations have 
found t hat the hot gas halos of sate l lites are not stripped in - 
stantly (|Kawata fc MulchaevjboO S': 'McCarthv et al.l l2008 



Recently, severa l authors (jKang fc van den Bosc 



I Font et al] |2008| ) have proposed improved recipes for the 
treatment of cooling onto satellites, which produce better re- 
sults for the predictions of satellite colours. Clearly, our anal- 
ysis should be repeated with one of these improved treat- 
ments implemented in our models. However, it is probably 
also important to properly treat the effects of ram pressure 
stripping on b oth the satellites' hot halo and the cold gas 
in the galaxv jLanzoni et al.l 20051: iQuilis. Moore, fc Bower! 



I2OOOI : lOkamoto fc NagashimalbOOa ). This will clearly be an 



important area for improvements to the modelling and fur- 
ther investigations. 
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APPENDIX : THE EFFECT OF DUST, 
SELECTION CRITERIA, AND GROUP 
CATALOG HALO MASS ESTIMATES 

All of our models include a treatment of dust extinction, 
which affects both the colours and magnitudes of the model 
galaxies. In addition, as mentioned in §4, for our main analy- 
sis we have applied selection criteria to the models to mimic 
those that we believe to be present in the SDSS group cat- 
alog. We include only galaxies with apparent r-band mag- 
nitude brighter than 17.77 mags, and we also excluded ha- 
los which do not contain any galaxy member brighter than 
'^'^Mr —19.5 -f 5 log ft. In this Appendix we provide the 
results for the 2-d distributions of /red and /passive for the 
theoretical models without dust extinction and without se- 
lection criteria applied. In addition, we test for possible bi- 
ases that may arise from the approach used to assign halo 
masses to the groups in the SDSS group catalogs by applying 
this method to the model mock catalogs. 

In Fig.|9]we show /red and /passive without any 'dust cor- 
rections applied to the models. Since we do not have magni- 
tude information without dust corrections for the dL06 mod- 
els, only the S08 fiducial and MORGANA models are shown. 
Interestingly, the results now look much more similar to the 
results seen before for /passive- In particular, the blue ridge 
which was visible in the /passive diagrams corresponded to 
a red ridge in the /rod diagram (Fig. [Sjl . The blue ridge is 
now visible in /red as well, indicating that these galaxies are 
actually actively star forming, and were predicted to be red 
only because of dust extinction. In the case of /passive, the 
dust correction only affects the galaxy selection and causes 
a negligible change in the diagrams. The treatment of dust 
extinction is one of the most uncertain aspects of the mod- 
elling, and this highlights the advantage of using intrinsic 
physical quantities extracted from the observations. 

Fig. [To] shows the model predictions with no selec- 
tion effects applied. Comparing this with Fig. |6l we see 
that the results appear unchanged above a halo mass 
logMhaio/ft"^M0 ^ 11.6 and a stellar mass Mgai ^ IO^Mq. 
This is reassuring in the context of our present analysis. 
However, we can also see that there is interesting predicted 
behavior at lower halo and galaxy masses than we can cur- 
rently probe, and also interesting differences between the 
models at these masses. This suggests that it would be ex- 
tremely useful to obtain similar data that are complete to 
fainter levels, so that we could probe lower mass galaxies 
and lower mass halos. 

Finally, we investigate the procedure used to assign 
halo masses to the groups that are identified in the SDSS 
group catalog. For the results present ed here, ste l lar m asses 
for several models are obtained from iBeU et al.1 l|2003h hke 
Y07, and halo masses are assigned based on the "charac- 
teristic" stellar mass of the group, where the characteristic 
stellar mass is defined as the total stellar mass contributed 
by galaxies with "'^Mr ^5 —19.5 + 5 log ft. Halo masses are 
then assigned by matching the rank-ordered list of group 
characteristic stellar masses with a rank-ordered list of dark 
matter halo masses from a theoretical estimate of the DM 
halo mass function, assuming a monotonic mapping between 
the characteristic stellar mass and DM halo mass (see >i2.1l 
and Y07) . We apply this procedure to the halos in the mock 
catalog produced with the S08 fiducial semi-analytic model. 
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Figure 9. The fraction of "red" (/rod) {^cft panels) and "passive" galaxies (/passive) (right panels) in models without dust eorreetions. 
The colour scale is as in Fig. \E\ We present the results for central galaxies (left) and satellite (right) galaxies seperately. We see that the 
details of the /^ed distribution are quite sensitive to the dust correction, whereas /passive is not noticably affected by the dust correction. 
We also note that the dust-free /^od results appear more similar to the /passive results, which presumably probe the physical properties 
of galaxies more directly. 
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Figure 10. The fraction of "red" galaxies (/red) (left panels) and of "passive" galaxies (/passive) (right panels) without selection criteria. 
The colour scale is as in Fig. |6] We present the results for central galaxies (left) and satellite (right) galaxies seperately. It can be inferred 
that our main results do not depend on our selection criteria. 



and show the comparison between the true and estimated 
halo mass in Fig. 1111 We see that although the mean halo 
mass is estimated fairly accurately, there is a large scatter in 
true halo mass at a given value of the estimated halo mass. 

In Fig. 1121 we show again the distribution of /rod with 
halo mass and stellar mass, now using the Y07-like halo mass 
estimates for the semi-analytic models instead of the true 
halo masses. We see that the procedure artificially reduces 
the scatter in Mgai at fixed halo mass. This can be attributed 
to two effects. First, the halo mass estimates in the Y07 
group catalogue are based on the total stellar mass in the 
halo, which strongly correlates with the stellar mass of the 



central galaxy. In addition, the Y07 group catalogue includes 
corrections for various incompleteness effects in the SDSS 
(the factor C in Eqs. 3-4 in Y07), which creates scatter 
that is not visible in the mock catalogues. It should also 
be noted that the halo-mass estimating algorithm of Y07 
washes out many of the detailed features that are visible in 
the model predictions. Hence we should not make too much 
of the detailed features in the empirical Mgai-Mhaio diagrams 
but should focus on the mean trends instead. 
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Figure 11. The true halo mass in the fiducial semi-analytic model 
of SOS, compared with the halo mass estimate based on the total 
stellar mass in the halo, using a procedure similar to that of Y07 
for the SDSS-based group catalog. The colour-scale indicates the 
stellar mass of the central galaxy in each halo, as shown by the 
key on the figure. Although the mean halo mass is reproduced 
fairly well, there is quite a large scatter in the true halo mass at 
a given estimated halo mass. 
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Figure 12. The fraction of "red" galaxies (/red) ^ ^ function of 
halo mass and stellar mass, where halo masses have been assigned 
in the semi-analytic models using an approach similar to that used 
in the SDSS group catalogs. The colour scale is as in Fig. |6] We 
present the results for central galaxies (left) and satellite (right) 
galaxies seperately. We see that the procedure used to assign halo 
masses in the SDSS group catalogs reduces the scatter in M^^i at 
fixed halo mass, and washes out many of the detailed features of 
the 2d distribution that are visible in the raw model predictions. 



